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It will be discussed the statistics of the extreme values in time series characterized by finite-term 
correlations with non-exponential decay. Precisely, it will be considered the results of numerical 
analyses concerning the return intervals of extreme values of the fluctuations of resistance and 
defect-fraction displayed by a resistor with granular structure in a nonequilibrium stationary state. 
The resistance and defect-fraction are calculated as a function of time by Monte Carlo simulations 
using a resistor network approach. It will be shown that when the auto-correlation function of 
the fluctuations displays a non-exponential and non-power-law decay, the distribution of the return 
intervals of extreme values is a stretched exponential, with exponent largely independent of the 
threshold. Recently, a stretched exponential distribution of the return intervals of extreme values 
has been identified in long-term correlated time series by Bunde et al. (2003) and Altmann and 
Kantz (2005). Thus, the present results show that the stretched exponential distribution of the 
return intervals is not an exclusive feature of long-term correlated time series. 
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I. INTRODUCTION 

The puzzling questionsposed by the observations of 
climate changes P, [H, H, HI IE H have recently pro- 
vided further motivations to the investi gati on of extreme 
events and dynamics [1, [t| Ell EH Ha Il3|. In particu- 
lar, the statistics of the return times (RTs) of extreme 
values is one of the subjects that have recently attracted 
the attentions of many authors 0, EE El [lj], El El- 
Return times (or intervals) r qi of the extreme values of a 
fluctuating quantity x(t), associated with the overcoming 
of a given threshold q, can be defined d EH as the time 
intervals between two consecutive occurrences of the con- 
dition x(t) > q . In the case of uncorrelated events it is 
easy to show that the RTs are exponentially distributed 
[1 El EH ■ Exponential distributions describe the statis- 
tics of return intervals of prices in financial markets, wind 
speed data or daily precipitations in a given place for the 
same time windows [1 EH, E3 ■ 

On the other hand, in recent years it has been real- 
ized that several other important examples of time series 
display long-term correlations [H flt | \y j\- This is the 
case of phy siological data (heartbeats [20|,[21| and neuron 
spikes 1221), hydro- meteorological records (daily temper- 
atures [II El E3, III]), geophysical or astrophysical data 
(occurrence of earthquakes [H Hf| or solar flares W 



internet traffic [15| and stock market volatility [17J, [27 1 
records. Long-term correlated time series are character- 
ized by an auto-correlation function, C x (s), decaying as 
a power-law: 
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where the correlation exponent j x ranges between and 
1 (normalized series {xi}i=i....7v with zero average and 
unit variance are considered). Therefore, in these condi- 



tions, the mean correlation time r diverges (r is defined 
as the integral over s of C x ). The effect of long-terms 
correlations on the RT statistics has been recently stud- 
ied by Bunde et al. [H EB| and by Altmann and Kantz 
[ItJ ■ Both these studies were performed on the ground of 
numerical calculations on Gaussian and long-term corre- 
lated time series generated by the Fourier transform tech- 
nique and by imposing a power-law decay of the power 
spectrum 28]. The main conclusions of Bunde at al. 
were the following, i) The mean return interval R q is 
unchanged by the presence of long-terms correlations, ii) 
The distribution of the RTs becomes stretched exponen- 
tial: 



P q {r)=a/R q exp[-(6 r/i? 9 ) 7 ] 



(2) 



where the exponents 7 and "f x were found to be equal, 
iii) The RTs themselves are long-term correlated with a 
correlation exponent 7' ~ j x . Actually, as noted by Alt- 
mann and Kantz [l7j], the statement i) is true in general 
and can be identified with Kac's lemma introduced 
in the context of dynamical systems while the statement 
iii) is true only in very special conditions fl7| . Moreover, 
it must be noted that the results ii) and iii) only apply 
to linear time series (i.e. to series whose properties are 
completely defined by the power spectrum and by the 
probability density) [H E3- However, apart from this 
restriction, the stretched exponential distribution of the 
RTs can be considered a general feature in p resence of 
long-term correlations in a time series [H E3 ■ This fea- 
ture has important consequences on the observation of 
extreme events: indeed it implies a strong enhancement 
of the probability of having return intervals well below 
and well above R q , in comparison with the occurrence 
of extreme events in an uncorrelated time series. Fur- 
thermore, it must be noted [17| that the distribution in 
Eq. @ only depends on the parameter 7, being a and b 
function of 7. 
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Here it will be considered the effect on the distribu- 
tion of the RTs of the presence of finite-term correlations 
with non-exponential decay, a situation which can occur 
in systems which are approaching criticality, where inter- 
mediate behaviors with non-exponential and non-power- 
law decay of correlations can emerge [H| . Thus, in the 
next section it will be analyzed the RT statistics of ex- 
treme values of the resistance ad defect-fraction fluctua- 
tions displayed by a resistor with granular structure in a 
nonequilibrium stationary state [301 ] . 



II. METHOD AND RESULTS 

The time series analyzed consist in the fluctuations 
of resistance and defect-fraction of a thin resistor with 
granular structure, biased by an external current / and 
in contact with a thermal bath at temperature T. The 
resistance values are calculated by using the multi-species 
resistor network (MSN) model [30[. The MSN model de- 
scribes a conducting film with granular structure as a 
resistor network made by different species of elementary 
resistors. Precisely, the resistor species are distinguished 
by their resistances and by their energies associated with 
thermally activated stochastic processes of breaking and 
recovery of the elementary resistors. The broken elemen- 
tary resistors, named defects, correspond to resistors of 
resistance about 10 9 higher than that of normal resistors. 
As a result of the competition between pairs of breaking 
and recovery processes for each species, the network can 
reach a stationary state whose properties will be depen- 
dent on the external conditions. The resistance and the 
defect fraction of the network are calculated by Monte 
Carlo simulations as a function of the time for differ- 
ent temperatures [3(|. All the details about the MSN 
model and its results can be found in Ref. [3(|. It should 
be noted, however, that this model provides a correlation 
time of the resistance fluctuations dependent on the tem- 
perature: t ~ (T — T*)~ 9 with = 2.7, thus diverging 
for T — > T* . As a result, the power spectral density of 
the resistance fluctuations shows a 1//" dependence of 
frequency, with a ~ 1 for T < T* . It should be also 
noted that the MSN model represents an extension to 
systems characterized by 1 // noise of a previous model 
existing in the literature, the stationary and biased resis- 
tor network (SBRN) model [U HH, HH ■ In any case, it 
must be underlined that, apart from the specific system 
described by the MSN model, the method used here for 
generating the time series can be also viewed as a pure 
numerical algorithm for generating numerical series with 
different and tunable correlation properties. 

Long R(t) and p(t) time series (typically made of 
1 v2 X 10 6 records) have been generated and analyzed 
for different temperatures. Here R indicates the resis- 
tance of the network expressed in Ohm and p the defect- 
fraction, while the time t is expressed in iteration steps. 
The analysis has been performed by considering normal- 
ized series: x(t) = (R(t)— < R >)/<r, where < R > is the 
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FIG. 1: Normalized x(t) signal (zero average and unit vari- 
ance). Here the original records consist of resistance fluctua- 
tions calculated at 400 K. The time is expressed in simulation 
steps. 
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FIG. 2: Auto-correlation functions of the resistance fluctua- 
tions calculated at 300 K (up triangles), 400 K (down trian- 
gles) and 500 K (diamonds). The grey solid curves represent 
the best- fit with the function given in Eq. @, while the grey 
dashed curve shows the best-fit with a power-law. 



average value of the network resistance and a the root- 
mean-square deviation from the average (a similar pro- 
cedure has been adopted for the p(t) series) . The follow- 
ing quantities have been considered: the auto-correlation 
function, the PDF of the x records, the return intervals 
r q of the extreme values for different threshold q and 
their distribution P q (r) (q is expressed in units of a). All 
the results shown here are obtained by considering a net- 
work of size 75 x 75, biased by a current / = 1.0 A and 
made by 15 species of resistors. The details concerning 
the network can be found in Ref. [30l 

Figure 1 displays a typical trend of x(t). Only a small 
portion of the total number of records, N — 2 x 10 6 , is 
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FIG. 3: Auto-correlation functions of the defect-fraction fluc- 
tuations calculated at 300 K (up circles), 400 K (squares) and 
500 K (crosses). The grey solid curves represent the best-fit 
with the function given in Eq. ©. 



shown. In this case x(t) represents the normalized resis- 
tance fluctuations calculated at a temperature of 400 K. 
The probability density of several x(t) series has been also 
analyzed. It has been found that the PDF of the resis- 
tance fluctuations presents a significant non-Gaussianity 
that progressively increases at hig her temperatures (go- 
ing from 300 to 500 K) ^HHHf, while the PDF of the 
defect-fraction fluctuations is Gaussian at all the temper- 
atures. 

Figure 2 reports the auto-correlation functions, C x (t), 
of three resistance fluctuation time series obtained for 
three different temperatures. Precisely, the three black 
curves correspond to a temperature of 300, 400 and 500 
K respectively going from the upper to the lower one. 
The grey dashed curve shows the best-fit to the C x data 
at 300 K with a power-law of exponent 7 = 0.22, while 
the grey solid curves represent the best-fit to the data at 
400 and 500 K with the function: 



The fitting parameters are Co = 1.23, h = 
1.42 x 10 4 for the data at 400 K and C 



(3) 

0.30 and u = 
= 0.98, ti = 



0.36 and v! = 1.51 x 10 3 for the data at 500 K. Many 
other functions have been also considered for the best- 
fit of the C x data. However, it has been found that the 
function in Eq. ^ optimizes the best-fit procedure with 
the minimum numbers of fitting parameters. By using 
Eq. ([3]), it is possible to calculate the correlation times 
characterizing the resistance fluctuations: r = 1051 at 
400 K and r = 148 at 500 K (see Ref. 30 for further 
details). Figure 3 shows the auto-correlation functions 
of the defect-fraction times series obtained at 300, 400 
and 500 K (black curves, again going from the upper 
to the lower one) and the best-fit to the C x data with 
Eq. Q (grey solide curves). The fitting parameters are 
C = 1.12, h = 0.14 and u = 5.96 x 10 4 for the data 



FIG. 4: Distribution of the return intervals for values above 
the threshold q = 2.5 (dashed black curve at the border of 
the grey area) obtained from the resistance fluctuations cal- 
culated at 300 K (long-term correlated records). For com- 
parison, the distribution of the return intervals for the same 
threshold value and obtained after shuffling the data is shown 
by circles. 




FIG. 5: Double-logarithmic plot of the normalized probability 
density of the return intervals for different threshold values q. 
The x{i) records are the same of Fig. 4, i.e. they present 
long-term correlations. The slope of the dashed line is 0.3. 



at 300 K, C = 1.19, hi = 0.19 and u' = 8.17 x 10 3 
for the data at 400 K and C ' = 1.08, h" = 0.23 and 
u" = 1.29 x 10 3 for the data at 500 K. The correlation 
times are: r = 1.40 x 10 4 at 300 K, r = 1653 at 400 K 
and t = 306 at 500 K. 

Let us start the analysis of the return intervals by con- 
sidering first the RT distribution for extreme values of 
the resistance fluctuations at 300 K. As pointed out by 
Fig. 2, in this case C x scales as a power- law over more 
than four orders of magnitude, indicating the existence 
of long-term correlations. Figure 4 reports the distribu- 
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FIG. 6: Double-logarithmic plot of the normalized probability 
density of the return intervals for different threshold values q. 
In this case the data (resistance fluctuations at 400 K) are 
no more long-term correlated but exhibit intermediate-term 
correlations. The slope of the dashed line is 0.26. 



FIG. 7: Double-logarithmic plot of the normalized probability 
density of the return intervals for different threshold values q. 
Also in this case the data (defect-fraction fluctuations at 400 
K) are characterized by intermediate-term correlations. The 
slope of the dashed line is 0.23. 



tion of r q for q = 2.5 in semi- logarithmic plot. Precisely, 
a normalized representation has been adopted for conve- 
nience, by reporting the product of mean return inter- 
val R q for the probability density P q {r) as a function of 
the ratio r q /R q . For comparison, the distribution of r q 
for the same threshold, obtained after shuffling the x(t) 
records is also shown. Thus Fig. 4 highlights the typical 
effect of long-term correlations: strong enhancement of 
the probability of having return intervals well below and 
well above R q , in comparison with the occurrence of ex- 
treme values in uncorrelated time series, as described in 
Refs. [TH, [l|| [l?]. Actually, the distribution of r q is found 
to be a stretched exponential, in agreements with the re- 
sults of Refs. [H], [l|| [l?]. This is shown in Fig. 5, which 
reports a double-logarithmic plot of the normalized prob- 
ability density of the r q as a function of the ratio r q /R q 
and for different thresholds. In this representation, the 
slope of the dashed line directly provides the exponent 7 
in Eq. ©. As expected [H, El, Il3| , 7 = 0.3 « j x . 

Now, let us consider the RT distribution of extreme 
values for time series with finite correlation time and 
slow, non-exponential decay of correlations: a situation 
that can be called of " intermediate-term correlations" . 
As shown by Figs. 2 and 3, the auto-correlation func- 
tions of the resistance fluctuations at 400 and 500 K (i.e. 
at T > T*) and of the defect-fraction fluctuations at 
300 -f- 500 K are well described by Eq. ©. Thus, these 
series just exhibit this intermediate behavior of the time 
correlations. Figures 6 and 7 display double-logarithmic 
plots of the normalized probability density of the r q val- 
ues calculated for the fluctuations at 400 K respectively 
of the resistance and of the defect-fraction. Different 
thresholds q ranging from 2 to 3.5 are considered in both 
cases. The slope of the solid line is 0.26 and 0.23, respec- 
tively for the case of Figs. 6 and 7. Figures 6 and 7 show 



that the distribution of the return intervals of extreme 
values of these series is well described by a stretched ex- 
ponential and that the value of the exponent 7 is inde- 
pendent of the threshold q in a large range of q-values. 
This occurs even in absence of long-term correlations and 
in presence of a finite correlation time. 



III. CONCLUSIONS 

The distribution of return intervals of extreme val- 
ues has been studied in several time series with different 
correlation properties: long-term and finite-term correla- 
tions. Precisely, it has been analyzed the return interval 
distribution of the fluctuations of resistance and defect- 
fraction displayed by a resistor with granular structure 
in noncquilibrium stationary states at different temper- 
atures T. The resistance fluctuations were calculated by 
using the MSN model which is based on a resistor net- 
work approach [30]. It has been found that when the 
auto-correlation function displays a non-exponential and 
a non-power-law decay, the distribution of the return in- 
tervals is well described by a stretched exponential with 
exponent 7 largely independent of the threshold q. This 
result shows that the stretched exponential distribution 
describes the distribution of the return intervals of ex- 
treme values not only when long-term correlations are 
present in the time series [lj| [l|| [Vq , but even when 
finite-term correlations with non-exponential decay exist 
among the records, a situation typical of many systems 
which are approaching criticality. 
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